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ABSTRACT 



We use the iBusca et al.l (120121 ) measurement of the Hubble parameter at 
redshift z = 2.3 in conj unction with 21 lowe r z me asure ments. fromlSimon et al. 



fl2005[ ). lGaztanaga et al.l fl2009[ ) . IStern et al.l fl2010f ). and lMoresco et all (120121 ). to 
place constraints on model parameters of const ant and time-evolvi ng dark energy 
cosmological models. The inclusion of the new iBusca et al.l ( 12012| ) measurement 
results in H(z) con straints significantly more restrictive than those derived by 
Farooq et al.l (|2012l ). These H[z) constraints are now more restrictive than those 



that follow from current Typ e la supernova (SNIa) apparent magnitude mea- 
surements (jSuzuki et al.ll2012l ). The H{z) constraints by themselves require an 
accelerating cosmological expansion at about 2 a confidence level, depending on 
cosmological model and Hubble constant prior used in the analysis. A joint anal- 
ysis of H{z), baryon acoustic oscillation peak length scale, and SNIa data favors 
a spatially-flat cosmological model currently dominated by a time-independent 
cosmological constant but does not exclude slowly-evolving dark energy density. 



1. Introduction 

Observations indicate the cosmological expansion is accelerating now and that the Uni- 
verse is spatially fiat, provided the dark energy density res ponsible f o r the a cceleration is clos e 
to or time in d epend e nt. For revie ws of dark energy see L 
Bolotin et al.l (12011! ) , IWangI (120121 ) , and references therein 



ET 



imened (I2OIII ). iLi et al.l (l2011al ). 



^ Instead of dark energy that dominates the current cosmological energy budget, a less likely possibility 
is that these observations are an indication that general relativity needs to be modified on cosmological 
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In the "standard" spatially- flat ACDM cosmological model ( lPeebleslll984j ) dark energy 
— Einstein's cosmological constant A — contributes a little more than 70% of the current 
energy budget. Non-relativistic cold dark matter (CDM) is the next largest contributer 
(at a little more t han 2 0%), followed by non-relativistic baryonic matter (around 5%). See 
Ratra fc Vogeleyl ( l2008l ) and references therein for reviews of the standard model. It has 
been known for a while that the sta ndard ACDM model is reasonably compatible with most 



observational con straints (see, e.g., IJassal et al.ll2010l : IWilson et al.ll2006l : iDavis et al.l 12007 



Allen et all l2008l )l^ In the ACDM model the dark energy density is constant in time and 
does not vary in space. 

It is well known that the standard ACDM model has some puzzling features that are 



time ( 


Peebles & Ratra 


19881: Ratra & Peebles 


19881). For recent discussions of time-varvine 


dark energy models see 


Shevkhi et al. 


(2OI2I). 1 


3rax & Davis 


(2012). HoUensi 


tein et al. (2012). 


Cai et 


al. (2012) 


Ponce de Leon (12012 


). Costa et al. (2012). 


Gu et al. (2012 


). Basilakos et al. 


(2012) 


. Xu et al. 


(2012 


1, and references therein. In this paper we study two dark energy mod- 



els (with dark energy being either a cosmological constant or a slowly-evolving scalar field 
(h) as well dark energy parameterization. 

In the ACDM model, time-independent dark energy density (the cosmological constant 
A) is modeled as a spatially homogeneous fluid with equation of state pa = — Pa that relates 
the fluid pressure and energy density. The XCDM parameterization has often been used 
to describe slowly- decreasing dark energy density. In this case dark energy is modeled 
as a spatially homogeneous X-fiuid with equation of state px = wxPx- Here px and px 
are the pressure and energy density of the X-fiuid and the equation of state parameter 
Wx < —1/3 is independent of time. When wx = — 1 the XCDM parameterization reduces 
to the complete, consistent ACDM model. For all other values of wx < —1/3 the XCDM 
parameterization is incomple te as it does not describe spatial inhomogeneities (see, e.g. iRatra 
I991I : iPodariu fc Ratrall2000l ). For computational simplicity, in the XCDM case we consider 
only a spatially-fiat cosmological model. 



The 0CDM model is the simplest, complete and consistent model of slowly-decreasing 



length scales. See lBolotin et al.l (|201l[ ). ICapozziello fc De LaurentisI (|201ll ). IStarkmanl (|201l[ ). and references 
therein, for reviews of modified gravity. Here we assume that general relativity is an adequate description 
of cosmological gravitation. 

^ Note, however, that there are tentative observational indications that the "standard " CDM structure 
formation model, which is a ssumed in the ACDM model, might need modification (e.g. 
2003t IPerivolaropoulosllioiol ). 



Peebles fc Ratra 
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dark energy density ([Peebles fc Ratral Il988l : iRatra fc Peebles! Il988l ) . Here dark energy is 
modeled as the gradually decreasing (in (p and time) potential energy density V{(f)) of the 
scalar field. In this paper we assume an i nverse power-law poten tial energy density ^(0) oc 
</)~", where a is a nonnegative constant ( iPeebles &: Ratral llQSSf ). When a = the 0CDM 
model reduces to the corresponding ACDM model. For computational simplicity, we assume 
a spatially-fiat cosmology for 0CDM. 



(Suzuki et al. 


2012; 


Li et al. 


2011b: 


Astier & Pain 



2012l : lRuiz etaDl2012 . and references therein) ; cosmic microway e background (CMB) anisotropy 



measurements (e.g.. |Podariu et al.ll2001bt iKomatsu et al.ll201l[) combined with low estimates 



of the cosmological mass density Jchen fc Ratral boosi and references therein), provided 



the dark energy density is close t o or time independent; and, baryon acoustic oscillation 



(BAO) peak length sc ale data (e.g.. iBeutler et al.ll201ll ; iBlake et al.ll201ll ; Mehta et al.ll2012 



Anderson et al.ll2012l ). Current error bars associated with these three types of data are still 
too large to allow for a significant observational discrimination between the ACDM model 
and the two simple time-varying dark energy models discussed above. Additional data are 
needed for this task, as well as to provide a cross check on the above results. 

Other data that have been used for this purpose include lookback ti r ae as a func- 
tion of redshift ( jSamushia et al.l |2010| ; iDantas et al.l 1201 It iTonoiu et al.l 1201 It iThakur et al. 
2OI2I . and references therein) , gamma-ray burst luminosity distance a s a function of redshift 



Samushia fc Ratra 



2OIOI 



f 



(e.g., 

gravitational lensing ( IChae et al. 



Wang &: Dai 



2004 



201 



Lee fc Ng 



Busti et al.ll2012[ |Poitr3|2012|), strong 



20071 ; IZhang fc Wull2010l ; ICao et al.l 12012 



and r eferences therein), HII starburst galaxy appareri t magnitude as a function of redshift 



(e.g., 
shift ( 


Plionis et al. 


2010, 


2011 




Mania & Ratra 


2012), angular 


size as a function of red- 
id references therein), and 


Chen & RatrE 


I 2012; 


Lima &; Cunha 


2012; 


Jackson 


2012. ai 


galaxy cluster properties ( 


e.g.. 


Campanelli et al. 


2012; 


Devi et al. 


2011; 


De Boni et al. 


2012; 



Gonzalez et al .1 120121 ). The constraints from these data are, at present, significantly weaker 
than those from SNIa, BAO, and CMB anisotropy measurements, but it is anticipated that 
future data of these kinds will provide significant constraintsjfl 

Two other current data sets provide interesting constraints on cosmological parame- 
ters, somewhat comparable to those from SNIa, BAO, and CMB anisotropy data. These 
are galaxy cluster gas mass fraction as a function of redshift measurements (e.g., I Allen et al. 



In addition to soon to be available CM B anisotropy data from Planck, future space-based SNIa, BAO 



like, and galaxy cluster measurements (e.g., 



Podariu et al 



_2001al: ISamushia et al.lboill: ISartoris et al-lboij 



Basse et al.ll2012t iMaierotto et al.ll2012t iPavloy et al.ll2012| ) should soon provide interesting constraints on 



cosmological parameters. 
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2008 




Samushia & Ratra 


2008: Tone k Noh 


2011 


: Lu et al. 2011: Solano & Nucamendi 2012) 


and measurements of the Hubble parameter as a function of redshift (e.g., 


Jimenez et al. 


2003 




Samushia & Ratra 


2006: Samushia et al. 


2007: 


Sen & Scherrei 2008: 


Chen & Ratra 


2011b 


Kumar 2012: Wane & ZliansJ2011: 


Duan et al. 


2011: Aviles et al. 2012: Seikel et al. 


2012h. Interestinelv. most measurements now provide largely compatible constraints on 



cosmological parameters that are consistent with a currently accelerating cosmological ex- 
pansion. This provides confidence that the broad outlines of a standard cosmological model 
are now in place. 



In this paper we use th e 21 H{z) measurements of 



Simon et al.l (120051 ). iGaztafiaga et al. 



( 2009 ). IStern et all (120101 ). and Moresco et al.l (120121 ) [listed i n Table 1 o f iFarooa et al. 
(120121 )]. in conjunction with the H{z = 2.3) measurement of iBusca et al.l (|2012| ). deter- 



mined from BAO in the Lya forest (in combination with WMAP CMB anisotropy data, 
Komatsu et al.|[2011 ). to const rain the ACDM and 0CDM models and the XCDM parametriza- 
tion. The inclusion of the newlBusca et al.l (120121 ) measurement results in tighter constraints 
than those derived by Farooq et al. ( 20121 ) from the 21 H{z) measurements alone. The new 
H(z) constraints derived here are more restrictive than those derived from the recent SNIa 



data compilation of ISuzuki et al.l ( 20121 ). which more carefully accounts for the systematic 
errors in SNIa measurementsQ In addition to deriving H{z)-data only constraints, we also 
use these H{z) measurements in combination with recent SNIa and BAO data to jointly 
constrain cosmological parameters. Adding the H{z) data tightens the constraints, quite 
significantly in some parts of parameter space. More precisely, the H{z) measurements more 
significantly tighten constraints on the nonrelativistic matter density parameter than on the 
parameter that more closely controls the time evolution of the dark energy density 

Our paper is organized as follows. In Sec. [2] we present the basic equations of the 
three dark energy models we consider. Constraints from the data are derived in Sec. O We 
conclude in Sec. jH 



2. Dark energy models 

In this section we list relevant characteristics of the two models (ACDM and 0CDM) 
and the one parametrization (XCDM) we use in our analyses of the data. 



The study of H{z) data is a much less-developed field than that of SNIa data, so it is not impossible 
that future H(z) error bars might be larger than what we have used in our analysis here. 
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In the ACDM model with spatial curvature the Hubble parameter evolves as 

H{z- Ho, p) = Ho [Qmoil + zf + ^]A + (1 - Q^o - ^A)il + z)'] , (1) 

where Ho is the current value of Hubble parameter (the Hubble constant), the current value 
of the spatial curvature density parameter is Qko = 1 — ^mo ~ ^a, and the model parameter 
set we want to constrain is p = (f2mo,^A)- Here Qmo is the nonrelativistic (baryonic and 
cold dark) matter density parameter and fl^ is the time-independent cosmological constant 
density parameter. Below we shall have need for the dimensionless Hubble parameter E[z) = 
H{z)/Ho. 

The XCDM parameterization Friedmann equation is 

Hiz; Ho, p) = i/o[^^™o(l + + (1 - ^^^o)(l + ^)3(i+-x)]i/2^ (2) 

where for computational simplicity we consider only flat spatial hypersurfaces and the model 
parameters p = [Qmo^^x)- The XCDM parametrization is incomplete, as it cannot describe 
the evolution of energy density inhomogeneities. 



The 0CDM model ( jPeebles fc Ratral Il988l ) is the simplest, complete and consistent 
dynamical dark energy model. In this model dark energy is a slowly-rolling scalar field 
with an, e.g., inverse-power-law potential energy density V{(f)) = Km^cf)'" where nip = 
is the Planck mass, G is the Newtonian gravitational constant, and a is a non- negative free 
parameter that determines k. The scalar field part of the (pCDM model action is 

ml f , /I 



S=T7^ I V^in ^""9,(1)3,(1) - Kml(f)~'' ] d^'x, (3) 



IQn J " " V2 
where g'^'^ is the metric tensor and a and k are related as 

a + 4\ /2a(a + 2)^°/^ 



with corresponding scalar field equation of motion 



4> + 3-4>- Kaml<j)-^'^+^^ = 0, (5) 
a ^ 

where the overdot denotes a derivative with respect to time. In the spatially-fiat case the 
Friedmann equation for the 0CDM model is 

H{z- Ho, p) = Ho[n„,o{l + zf + %{z, a)Y/\ (6) 
where 0.^{z, a) is determined by the field energy density 

"^p / 1 i 2 I 2 



167r V2 



+ Km^r'' . (7) 
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Fig. 1. — Thick solid (t hin dot-dashed) line s correspond to 1, 2, and 3 a constraint con- 
tours from the new (old, iFarooq et al.ll2012l ) H{z) data for the ACDM model. The filled 
(empty) circle is the best fit point from the new (old) H[z) data. The left panel is for 
the i^o = 68 lb 2.8 km s~^ Mpc~^ prior and the right panel is for the Hq = 73.8 ± 2.4 
km s~^ Mpc~^ one. The dashed diagonal lines correspond to spatially-flat models, the 
dotted lines demarcate zero-acceleration models, and the shaded area in the upper left- 
hand corners are the region for which there is no big bang. The filled circles correspond 
to best-fit pair {Qmo,^A) = (0.21,0.53) with Xmm — ^^.l (left panel) and best-fit pair 
(f2mo,^A) = (0.26,0.77) with xLin = 16.1 (right panel). The empty circles correspond 
to best-fit pair (f2mO;^A) = (0.28,0.62) with xLin = 14.6 (left panel) and best-fit pair 
(^]„o,fiA) = (0.42,0.97) with xLn = 14.6 (right panel). 



Equations ([5]) — (I7j) constitute a system of differential equations which can be solved numer- 
ic ally for the (/)CDM mod el Hubble parameter H{z), using the initial conditions described 



m 



Peebles fc Ratral (119881 ). In this case the model parameter set is p = {flmo, ct) 



3. Constraints from the data 



We first study H{z) data constrain ts on cosmological pa rameters. To the 21 indep e ndent 
H{z) data points listed in Table 1 of iFarooq et al.l (120121 ) we add the iBusca et al.l (120121 ) 
H{z = 2.3) = 224 ± 8 km s~^ Mpc~^ measurement, determined frorn BAO in the Lya forest 
(in conjunction with WMAP CMB anisotropy data. iKomatsu et al.ll201ll ). 
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Fig. 2. — Thick so lid (thin dot-dashed ) hnes correspond to 1, 2, and 3 a constraint contours 
from the new (old. iFarooq et aLll2012l ) H{z) data for the XCDM model. The filled (empty) 
circle is the best fit point from the new (old) H{z) data. The left panel is for the Hq = 68±2.8 
km s~^ Mpc~^ prior and the right panel is for the Hq = 73.8 ± 2.4 km s~^ Mpc~^ one. 
The dashed horizontal lines at ux = — 1 correspond to spatially- flat ACDM models and 
the curved dotted lines demarcate zero-acceleration models. The filled circles correspond 
to best-fit pair {flrno,'-^x) = (0.27,-0.82) with xLin = 15.2 (left panel) and best-fit pair 
(Qmo,oJx) = (0.36,-1.1) with xLm — 15-9 (right panel). The empty circles correspond 
to best-fit pair (Qrno^^x) = (0.31,-0.94) with xLin — 1^.6 (left panel) and best-fit pair 
{nmo,uJx) = (0.30, -1.30) with xLn = 14.6 (right panel). 



To constrain cosna ological parameters p of the models of interest we follow the procedure 
of iFarooq et al.l (|2012| ) . We again marginalize over the nuisance parameter Hq using two 
different Gaussian priors with Hq ± (Jho= 68 ± 2.8 km s~^ Mpc~^ and with Hq ± cthq^ '''3-8 
± 2.4 km s~^ Mpc~^. As discussed there, the Hubble constant measurement uncertainty 
can significantly affec t cosmological parameter estimation (for a recent example see, e.g.. 



Calabrese et al.l 120121). T he lower of the two values we use is from a me dian statistics 



analysis (iGott et al.ll200ll ) of 553 meas urements of Hq (Chen fc Ratrall2011al) : this estimate 



has been stable for over a d ecade now (IGott et al.ll2001 



Chen et al.ll2003l ). The other value 



is a recent, HST based one (iRiess et al.ll201ll ). Ot her recent estimates a r e compatible with 



at least one of the two values we use (see, e.g., iFreedman et al.l l2012l : ISorce et al.l 12012 



CoUess et al.ll2012l : iTammann fc Reindlll2012l ). 



We maximize the likelihood jCh(p) with respect to the parameters p to find the best- 
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Fig. 3. — Thick so lid (thin dot-dashed ) hnes correspond to 1, 2, and 3 a constraint contours 
from the new (old. iFarooq et al.ll2012l ) H[z) data for the 0CDM model. The filled (empty) 
circle is the best fit point from the new (old) H{z) data. The left panel is for the Hq = 68±2.8 
km s~^ Mpc~^ prior and the right panel is for the Hq = 73.8 ± 2.4 km s~^ Mpc~^ one. 
The horizontal axes at a = correspond to spatially-flat ACDM models and the curved 
dotted lines demarcate zero-acceleration models. The filled circles correspond to best-fit 
pair {Qrno,Ci) = (0.36,0.70) with xLin = 15.2 (left panel) and best-fit pair {Qmo,Ci) = 
(0.25,0) with Xmin = 16.1 (right panel). The empty circles correspond to best-fit pair 
(n^o,a) = (0.30,0.25) with xLn = 14.6 (left panel) and best-fit pair (fi„^o,a) = (0.27,0) 
with Xmin = 15.6 (right panel). 



fit parameter values po. In the models we consider = — 21n[£j:/(p)] depends on two 
parameters. We define 1, 2, and 3 a confidence intervals as two-dimensional parameter sets 



bounded by xlip) = xUpo) + 2.3, Xnip) = xUpo) + 6.17, and x%ip) 
respectively. 



xUpo) + 11.8, 



Figur es [H— [SI show the co nstraints from the H{z) data derived here, as well as those 
derived by lFarooq et al. J2012h . for the three da rk energy models, a nd for the two different Hq 
priors. Clearly, the H{z = 2.3) measurement of iBusca et al. J2012h significantly tightens the 
constrains. Given that the nonrelativistic matter density is larger relati v e to t he dark energy 
density at ;z = 2.3, it is perhaps not unexpected that the iBusca et al.l ( 120121 ) measurement 
tightens the constraints on Qmo much more significantly than it does for the constraints on 
the other parameter which more strongly affects the evolution of the dark energy density. 
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Fig. 4. — Thick solid (thin dot-dashed) hnes are 1, 2, and 3 a constraint contours for the 
ACDM model from a joint analysis of the BAO and SNIa (with systematic errors) data, with 
(without) the H{z) data. The full (empty) circle marks the best-fit point determined from the 
joint analysis with (without) the H{z) data. The dotted sloping line corresponds to spatially- 
flat ACDM models. In the left panel we use the Hq = 68 ± 2.8 km s~^ Mpc~^ prior. Here the 
empty circle [no H{z) data] corresponds to best-fit pair (i^mO; ^a) = (0.30, 0.73) with — 
551 while the full circle [with H{z) data] indicates best-fit pair (fimo^ ^a) = (0.29, 0.72) with 
Xmin — 567. In the right panel we use the Hq = 73.8 ± 2.4 km s~^ Mpc~^ prior. Here the 
empty circle [no H{z) data] corresponds to best-fit pair {Qmo, ^a) = (0.30, 0.73) with Xmin = 
551 while the full circle [with H{z) data] demarcates best-fit pair (fi^c^A) = (0.28,0.78) 

with Xmin = 568. 



see Figs. [2] and |3l 

Comparing the H{z) constraint s derived here , and shown in Figs. [U— [3] here, to the 
SNIa constraints shown in Fig. 4 of iFarooq et al.l (120121 ). we see that the new H[z) data 
constraints are significantly more restrictive than those that follow on using the SNIa data. 
This is a remarkable result. Qualitatively, because of the dependence on the Hq prior and on 
the model used in the analysis. Figs. [1]— [3] show that the H{z) data alone require accelerated 
cosmological expansion at approximately the two standard deviation confidence level. 



While the H{z) data provide tight constraints on a linear combination of cosmological 
parameters, the banana-like constraint contours of Figs. [I]— [3] imply that these data alone 
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Fig. 5. — Thick solid (thin dot-dashed) hnes are 1, 2, and 3 a constraint contours for 
the XCDM parametrization from a joint analysis of the BAO and SNIa (with systematic 
errors) data, with (without) the H{z) data. The full (empty) circle marks the best-fit point 
determined from the joint analysis with (without) the H{z) data. The dotted horizontal line 
at cux = — 1 corresponds to spatially- flat ACDM models. In the left panel we use the Hq = 68 
± 2.8 km s"^ Mpc-i prior. Here the empty circle [no H{z) data] corresponds to best-fit pair 
{flmo,uJx) = (0.30, —1.03) with Xmin = 551, while the full circle [with H{z) data] demarcates 
best-fit pair {Qrno,^x) = (0.29,-0.99) with Xmin = 568. In the right panel we use the Hq 
= 73.8 ± 2.4 km s~^ Mpc~^ prior. Here the empty circle [no H{z) data] corresponds to 
best-fit pair {Qjno,'^x) = (0.30, —1.03) with xLin — 551 while the full circle [with H{z) data] 
indicates best-fit pair (fi^ci^x) = (0.28, —1.05) with xLin — 569. 



cannot significantly discriminate between cosm ological model s . To tighten the constraints 
we must add other data to the mix. Following iFarooq et al.l (120121 ). we derive constraints 
on cosmological parameters of the three models from a joint analysis of the H(z) da t a wit h 
the 6 BAO peak length scale measurements of iPercival et al.l (120101 ) , iBeutler et al.l (120111 ) , 
and iBlake et al.l (l201ll ). and the Union2.1 compilation of 580 SNIa apparent na agnitude 



measurements (covering a redshift range 0.015 < z < lA) from lSuzuki et al.l (|2012| ). 

Figures imn] show the constraints on cosmological parameters for the ACDM and 0CDM 
models and the XCDM parametrization, from a joint analysis of the BAO and SNIa data, 
as well as from a joint analysis of the BAO, SNIa, and H{z) data. Including the H{z) data 
in the analysis tightens the constraints, somewhat significantly (sometimes by more than 
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Fig. 6. — Thick solid (thin dot-dashed) hnes are 1, 2, and 3 o constraint contours for the 
0CDM model from a joint analysis of the BAO and SNIa (with systematic errors) data, 
with (without) the Hi^z) data. The full (empty) circle marks the best-fit point determined 
from the joint analysis with (without) the Hi^z) data. The a = horizontal axes correspond 
to spatially-flat ACDM models. In the left panel we use the ifo = 68 ± 2.8 km s~^ Mpc~^ 
prior. Here the empty circle corresponds to best-fit pair (fimO; Q^) = (0.30, 0) with x^^^^ = 551 
while the full circle indicates best-fit pair (fimO; ct) = (0.29, 0) with x^^j^ = 567. In the right 
panel we use the i?o = 73.8 ± 2.4 km s~^ Mpc~^ prior. Here the empty circle corresponds 
to best-fit pair (nmO;«) = (0.30,0) with x^jj^ = 551 while the full circle demarcates best-fit 
pair (fi^o, «) = (0.27, 0) with Xmin = 569. 

two standard deviations), in parts of the parameter spaces. Figure [7] shows the Hi^z) data 
and the two best-fit ACDM models. The ii{z^ data do support the idea of a deceleration to 
acceleration transition somewhere in the range 0.5 < z <\. 

Table 1 lists the two standard deviation bounds on the individual cosmological param- 
eters, determined from their one- dimensional posterior probability distributions functions 
(which are derived by marginalizing the two-dimensional likelihood over the other cosmolog- 
ical parameter). 
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Fig. 7. — Measurements and predictions for H{z)/{1 + 2;) as a function of z. Dashed 
(dotted) lines show the predictions for the best-fit ACDM model from the combined 
BAO, SNIa, and H{z) data analyses, with cosmological parameter values {flmo, ^a, h) = 
(0.29, 0.72, 0.68)[(0.28, 0.78, 0.738)]. 



Conclusion 



Adding the iBusca et al.l ( 12012| ) z = 2.3 measurement of t he Hubble paramet er, from 
BAO in the Lya forest, to the 21 H{z) data points tabulated in iFarooq et al.l (|2012[ ). results 
in an H{z) data set that provides quite restrictive constraints on cosmologica l parameters 
These constraints are tighter than those that follow from the SNIa data of ISuzuki et al. 
( 120121 ). which carefully accounts for all known systematic uncertainties. The H[z) field is 
much less mature than the SNIa one, and there might be some as yet undetected H{z) 
systematic errors that could broaden the H{z) error bars, as has happened in the SNIa case. 
However, we emphasize that the observers have done a careful analysis and the error bars 
we have used in our analysis have been carefully estimated. In addition to providing more 
restrictive constraints, the H{z) data alone requires accelerated cosmological expansion at 
the current epoch at approximately 2 a confidence level, depending on model and Hq prior 
used in the analysis. 
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Model and prior 


BAO±SNIa 


BAO±SNIa±iJ(z) 


ACDM, h = 


- 0.68 


±0.028 


0.25 < QrnO < 0.36 

0.53 < i^A < 0.89 


0.26 < QrnO < 0.33 

0.60 < iZA < 0.84 


ACDM, h = 


0.73^ 


i ± 0.024 


0.25 < Qmo < 0.36 
0.53 < i^A < 0.89 


0.25 < ^]™o < 0.32 
0.66 < ilA < 0.89 


~\7' — \ TV T 7 

XCDM, h = 


-- 0.68 


± 0.028 


0.30 < QrrM < 0.38 
— i.io ^ u7x ^ — U. ( o 


0.27 < QrnO < 0.32 
1 IT? ^ / 1 ^ n 77 


XCDM, h = 


0.73 J 


i ± 0.024 


0.30 < Qmo < 0.38 
-1.18 < wx < -0.78 


0.25 < QrnO < 0.30 

-1.15 <ujx< -0.90 


0CDM, h = 


: 0.68 


±0.028 


0.25 < Qmo < 0.35 
< a < 0.54 


0.25 < Qrr^O < 0.32 

< a < 0.56 


0CDM, h = 


0.73^ 


i ± 0.024 


0.25 < Q^o < 0.35 
< a < 0.54 


0.25 < Qmo < 0.30 
< a < 0.21 



Table 1: Two standard deviation bounds on cosmological parameters using BAO±SNIa and 
BAO±SNIa±iJ(z) data, for three models and two Hq priors. 

In summary, the results of the joint analysis of the H{z), BAO, and SNIa data are quite 
consistent with the predictions of the standard spatially-flat ACDM cosmological model, with 
current energy budget dominated by a time-independent cosmological constant. However, 
currently-available data cannot rule out slowly-evolving dark energy density. We anticipate 
that soon to be available better quality data will more clearly discriminate between constant 
and slowly-evolving dark energy density. 

We thank Data Mania and Larry Weaver for useful discussions and helpful advice. 
This work was supported in part by DOE grant DEFG03-99EP41093 and NSF grant AST- 
1109275. 
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